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We study the time evolution of the distribution functions for hot and or degenerate gases of two 
flavors of Dirac neutrinos as a result of flavor mixing and dephasing. This is achieved by obtaining the 
time evolution of the flavor density matrix directly from quantum field theory at finite temperature 
and density. The time evolution features a rich hierarchy of scales which are widely separated in the 
nearly degenerate or relativistic cases and originate in interference phenomena between particle and 
antiparticle states. In the degenerate case the flavor asymmetry AN(t) relaxes to the asymptotic 
limit AN(oo) — AiV(O) cos 2 (2$) via dephasing resulting from the oscillations between flavor modes 
that are not Pauli blocked, with a power law 1/t for t > t s ~ 2Icf/AM 2 . if is the largest of 
the Fermi momenta. The distribution function for flavor neutrinos and antineutrinos as well as 
off-diagonal densities are obtained. Flavor particle-antiparticle pairs are produced by mixing and 
oscillations with typical momentum k ~ M the average mass of the neutrinos. An effective field 
theory description emerges on long time scales in which the Heisenberg operators obey a Bloch-type 
equation of motion valid in the relativistic and nearly degenerate cases. We find the non-equilibrium 
propagators and correlation functions in this effective theory and discuss its regime of validity as 
well as the potential corrections. 

PACS numbers: 14.60.Pq,12.15.Ff,11.90.+t 



I. INTRODUCTION 

Neutrinos are the bridge between particle physics, astrophysics, cosmology and nuclear physicspj E3- 0- E3l - and 
after almost four decades of the prescient suggestion that neutrinos may oscillatel3, [a|. a wealth of experimental data 
confirms that neutrinos are massive and that different flavors mix and oscillate flOl ITTl IT^L IT^| . Neutrino masses 
and mixing decidedly points to new physics beyond the standard model and profoundly impacts on the physics, 
astrophysics and cosmology of neutrinos. Neutrino oscillations in matter may provide an explanation of the solar 
neutrino problem by the resonant conversion of flavor neutrinos in the medium, namely the MSW effect [3 E3 
(for recent reviews see The dynamical aspects of neutrino oscillations in extreme conditions of temp~atu7e 

and density play an important role in Big Bang Nucleosynthesis (BBN) and in the lepton asymmetry in the early 
Universe 1 18j| (for a recent review seeQ) as well as in the physics of core collapse supernovae and the formation, 
evolution and cooling of neutron stars jllA l2fj l2l| . The study of the dynamical evolution of a hot and/or dense gas of 
neutrinos that include mixing as well as collisions has been and continues to be the subject of much attention in the 
literature. Neutrino mixing and oscillations introduce a novel aspect in the description of flavor equilibration, since 
the weak interactions involve flavor (weak) eigenstates while time evolution is described in terms of mass eigenstates. 
Therefore in a dense and/or hot medium where neutrino interactions cannot be neglected collisional processes must 
be studied on the same footing as the dynamics of oscillations. Furthermore in a dense background of neutrinos such 
as is the case in the early Universe during the time relevant for BBN or during the time scale of neutrino trapping in 
a protoneutron star, the neutral current interaction leads to a contribution to the neutrino self-energy from forward 
scattering off the neutrino background akin to the contribution from the electron plasma that leads to MSW resonance 
enhancement [T^| . In dense neutrino gases, this self-energy contribution leads to a non-linear problem for the evolution 
of a given neutrino interacting with the neutrino background. 

The dynamics of neutrino oscillations was originally studied in terms of Bloch-type equations akin to the equation 
of motion for a spin in a magnetic fieldP,|2,|3,0,|23| which are generally valid for single particle descriptions in the 
relativistic limit. For the case of single particle states this equation of motion for neutrino oscillations was derived 
from the underlying field theory in the relativistic limit |5l l23j. This formulation of the dynamics of oscillations of 
single particle states was extended to a kinetic description of oscillations and mixing in a medium pi, l25l l2f| . The 
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resulting equations in principle include the effects of collisions as well as the non-linearities arising from neutrino 
forward scattering off a neutrino background. They have been implemented to study the evolution of the neutrino 
distribution functions in the early Universe^, 113, HE |2j| HH El, m supernovaejsE HE E3 as weu as to study the 
relic neutrino asvmmetrv|35| . Novel fascinating self-synchronization phenomena emerges as a consequence of the 
non-linearities in a neutrino background with potential implications on CP (and baryon) asymmetry in the early 
Universe 0, HI] . 

An alternative quantum field theory treatment of neutrinos in the medium used the ingredients of thermal field 
theorv[36l l37| combined with a self-consistent treatment in the case of a neutrino background [3^. Since the main 
method in this approach relies on the equilibrium description of thermal field theory, there is an underlying assumption 
that the neutrino background is nearly in equilibrium. 

More recently the validity of the single particle picture that underlies the kinetic equations for neutrinos in a 
medium has been critically re-examined [39j. 

In our view, the study of neutrino oscillations and mixing in the case of a dense and or hot neutrino background either 
via the set of kinetic equations [3 HE HE| or the thermal field theory approach invoke a variety of approximations 
some of which are not very clear. In the kinetic description several approximations are involved, from neglecting 
interference terms between particles and antiparticles by restricting the Hamiltonian[24j to some time averaging and 
restriction to single particle evolution [2(j. Some of these approximations motivated the study of ref . |39| . 

A full quantum field theory treatment of neutrino mixing reveals a more complex picture of oscillations beyond that 
of the single particle description [40l Ell |42| . The authors of these references pointed out that a careful treatment of 
the Fock representation of flavor states leads to novel contributions to the oscillation formula even for single particle 
states. While it has been argued recently that Fock states of flavor neutrinos may not be relevant for S-matrix 
processes a quantum statistical mechanics of dense and or hot flavor neutrino gases must necessarily rely on the 
Fock representation (occupation number) for flavor neutrino states. 

A quantum statistical description of a dense and or hot gas of flavor neutrinos requires the notion of an occupation 
number which inevitably implies a description in terms of Fock flavor states. Furthermore a chemical potential 
associated with a flavor neutrino is a variable conjugate to the number of these flavor neutrinos. 

Regardless of whether the variety of approximations usually invoked are justified for practical purposes, the study 
of the dynamics of neutrino mixing and oscillations from the point of view of quantum field theory is clearly of 
fundamental importance as a prelude towards physics beyond the standard model. While there have been studies of 
the quantum field theory aspects in vacuum we are not aware of any previous study of the quantum field theory of 
mixing in a dense and or hot medium with neutrinos. 

The goal of this study: 

In the presence of flavor mixing, individual flavor number is not conserved and a density matrix that is diagonal in 
the flavor Fock basis will evolve in time and develop off diagonal elements. 

Hence time evolution of a dense or hot neutrino gas has to be studied as a quantum mechanical initial value problem: 
an initial density matrix which is diagonal in the flavor basis is evolved in time with the full Hamiltonian with flavor 
mixing. In this article we focus on studying precisely the time evolution of a dense or hot flavor neutrino gas in the 
simplest case of free field theory. Our goal is to study the evolution of an initially prepared density matrix which 
is diagonal in the flavor basis and describes a quantum gas of flavor neutrinos at finite density or finite temperature. 
We undertake the study of the dynamics in free field theory as a prelude towards a complete understanding of 
oscillation phenomena in weak interactions. The first step of any systematic program must be the understanding at 
the simplest level. As will be detailed below, studying the dynamics of oscillations and mixing in a dense and/or hot 
medium even at the level of free field theory reveals a wealth of subtle and important phenomena which leads to a 
firmer understanding of the validity of the various approximations as well as highlighting the potential corrections. 

The problem that we study can be stated succinctly as follows 1 : Consider that at a given initial time we have a 
"box" that contains a hot or dense gas of flavor neutrinos with a given single particle distribution consistent with 
Fermi-Dirac statistics , how does this ensemble evolve in time?, how do the populations of flavor neutrinos evolve in 
time?, how do flavor neutrinos propagate in the medium??. 

While our ultimate goal is to study the evolution in the presence of the weak interactions, we begin our study in 
this simplest free field theory case and the case of two flavors with the following goals in mind 

• To study the evolution directly from the underlying quantum field theory without making any approximations. 
This study will clarify the nature of the various approximations invoked in the literature and exhibit the potential 
corrections. 
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• By keeping the full evolution, the different time scales will emerge thus paving the way to providing a firmer 
understanding of coherence effects as well as the time averaging implied by several approximations. 

• A first principle derivation of kinetic equations and or Boltzmann equations require the propagators for the 
fields 44] in the medium. Thus the study of the evolution in free field theory is the starting point for a systematic 
treatment of oscillations and collisions in a medium with a neutrino background. 

• As it will become clear below, the study of even the simple free field theory case reveals a wealth of phenomena 
as a consequence of flavor mixing, which to the best of our knowledge has not been recognized and explored fully 
before in the case of finite temperature and density. The full quantum field theory treatment unambiguously 
reveals all the complexities associated with flavor mixing and allows a systematic implementation of several 
approximations which clarify the regime of validity of the single particle description and provide an understanding 
of the corrections. 

Brief summary of the results: Our main results are briefly summarized as follows, 

• The dynamics of neutrino oscillations of a dense and or hot gas of flavor neutrinos features a hierarchy of time 
scales. The fast time scales are associated with interference effects between particle and antiparticle states while 
the slow scales emerge from interference between particle states (or antiparticle states) of different masses. In 
the nearly degenerate or relativistic case the scales are widely separated and processes which involve interference 
between particle and antiparticle states become subdominant on the slow dynamics. 

• An initial flavor asymmetry relaxes towards an asymptotic value AiV(oo) = AiV(O) cos 2 (20) (with 6 being 
the mixing angle) with a power law oc 1/i as a consequence of dephasing. Pauli blocking manifests in that 
neutrinos of one flavor can only oscillate into un-occupicd states of neutrinos of different flavor and dephasing is 
a consequence of oscillations between Pauli unblocked flavor states. We obtain the explicit time evolution of the 
distribution functions as well as off-diagonal correlation functions. We discuss the phenomenon of flavor pair 
production by mixing and oscillations. This is a consequence of the overlap between particle and antiparticle 
states and results in the production of pairs of flavored neutrinos with typical momenta k ~ M, the average 
mass of the neutrinos. 

• In the nearly degenerate case (as suggested by the recent combined observations) or in the relativistic case as 
is likely to prevail in the early universe as well as in core collapse supernovae, the different time scales are 
widely separated. This allows to establish an "effective" (free) field theory description valid on the slow time 
scales. The equations of motion for Heisenberg operators in this effective description are the oft quoted Bloch- 
type equations, but the effective field theory also describes the quantum fields. This effective theory allows to 
construct the Feynman propagators which feature distinct non- equilibrium aspects and to clearly identify the 
potential corrections and is valid both in the relativistic as well as in the nearly degenerate case. 

Our study is organized as follows: in section II the theory corresponding to two flavors of neutrinos as well as the 
density matrix that describes an initial state of flavor neutrinos is presented. In this section we address the quantization 
aspects and point out the source of subtle mixing phenomena between particles and antiparticles, confirming previous 
results in the liter ature^fj- I n sections III and IV we study the evolution of the flavor asymmetry as well as that of 
the individual distribution functions focusing on the emergence of a hierarchy of scales and extracting the asymptotic 
long time dynamics as well as the phenomenon of flavor pair production via oscillations. In section V we present the 
"effective" field theory that describes the long-time dynamics and discuss its regime of validity. In this section we 
obtain the Feynman propagators and discuss their non-equilibrium aspects. In section VI we discuss the regime of 
validity of the several approximations as well as caveats in the formulation and present our conclusions. 

II. NEUTRINO MIXING AND FLAVOR DENSITY MATRIX 

We focus our attention on the evolution of Dirac neutrinos postponing the case of Majorana neutrinos for further 
discussion elsewhere. Furthermore, we restrict the discussion to the case of two flavors which provides the simplest 
scenario. Most of the results can be extrapolated to the case of three active flavors including the case of sterile 
neutrinos, but for the subtleties associated with CP violating phases which of course are of great interest but will not 
be addressed here. We will call the flavors the electron and muon neutrino, but the results apply more broadly to 
active-sterile oscillations. 
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Consider the Dirac neutrino fields with the Lagrangian density given by 



L = v e (x)(i0)v e (x) + v^xXiflM*) + ( *M ) ( ^ ) ( ffij) 



(11.1) 



where m e „ is the mixing and we have absorbed a potential phase into a field redefinition. The mass matrix can be 
diagonalized by introducing a rotation matrix such that 



u e (x) \_(C S\ ( Vi(a:) 



(11.2) 



where for simplicity of notation we defined 



C = cos(9; 5 = sin (9 (II.3) 
where 9 is the mixing angle. The diagonalized mass matrix then reads 



Mi \ _ ( C -S \ m e m efl \( C S 
M 2 J [ S C \ m e a m„ -S C 



(II.4) 



In the mass eigenstate basis, the Lagrangian density becomes 



C = &(a;)(i0 - Mi)^i(a:) + ^(as)(*^ - M 2 )^ 2 (x). (II.5) 

In what follows, we reserve the latin label i = 1,2 for the fields associated with the mass eigenstates i/j and the 
greek label a = e, /x for the fields associated with the flavor eigenstates v. 

Upon quantization in a volume V, the flavor field operators v a (x) at time t = are written as 



A 

where the index A refers to the Dirac spin index and we have kept the same notation for the field and its spatial 
Fourier transform to avoid cluttering of notation. A flavor Fock representation is defined by choosing the spinors U 
and V respectively. In principle these spinors can be chosen to be the positive and negative energy solutions of a 
Dirac equation with an arbitrary mass, in what follows we will choose these to be m e ; rriu, namely the masses of the 
flavor eigenstates in the absence of mixing. While we consider this to be a physically motivated choice, it is by no 
means unique and different alternatives have been discussed in the literature [40t l4ll l42j . 
Thus the spinors U and V are chosen to be solutions of the following Dirac equations 



7° (7 ■ k + m a ) = -u a {k) (II.7) 



Wa(*) = \A 2 + ml (11.8) 
The Dirac spinors U and V, are normalized as follows (no sum over the index a) 

and the creation and annihilation operators , ; /?g , obey the usual canonical anticommutation relations. 

On the other hand, upon quantization the field operators ipi( x ) associated with mass eigenstates at time t = are 
given by 
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A 

where the spinors -F, G are now solutions of the following Dirac equations 



7 ° (7 • fc + Mi) A = -^(fc) G w - A (11.11) 



fc,A 

Ei(k) = yjk 2 + M? (11.12) 
with the normalization conditions (no sum over the label i) 



n(*)t_p(i) _ (jWtfjW _ J . ^)tgf(i) 
A:, A fc,A' fc.A fc,A' 1 1 fc.A —A 

Similarly, the operators a and 6 satisfy usual canonical anticommutation relations 



A. Hamiltonian and Charges 

The total free field Hamiltonian for mixed neutrinos in the diagonal (mass) basis is given by 



H = Y / [Mk)(i -k + MMk)] = £ (afx a Z + b ¥J b nl J ) El{k) ' (IL14) 

k,i k,\,i 

Therefore the time evolution of the operators a, b is given by 

fc,A v ; fe,A 

bf (t) = bf e-^W*. (11.15) 

/c , A , A 

The free field Lagrangian density is invariant under independent phase transformations of the fields ^1,2, 

hence the individual U(l) charges 

Qi= J d 3 x tm*, *)=E RMa - ^Mi + x ] ( IL16 ) 

fe,A 

are time independent. 

The discussion that follows will focus on describing a statistical density matrix which is diagonal in the flavor basis 
and describes a hot and or dense ensemble of flavor neutrinos. This discussion requires the flavor Hamiltonian which 
is obtained from the Lagrangian density Ull.lfl for vanishing mixing m e a — 0, namely 

H, =H e + H, = J2 [M*)(7 ' * + m a )u a (k)] = VJ + - l) Wa (A), (11.17) 

k,a fc,A,a 

The flavor Hamiltonian above is invariant under independent phase transformations of the flavor fields v a , thus the 
individual flavor charges commute with Hf 

qa = J vi{x)v a {x) = J2 vt(k)u a <k) = J2 [«gS + 4?A - 45 t/3 S + x ] (IL18) 
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Using the transformation law 1|II.2(1 between flavor and mass eigenstates it is straightforward to find that the total 
charges are the same, namely 

J2^l(k,t)^i(k,t) = ^4(fc,iK(fc,f) Qi + Q 2 = q e + q t , (H-19) 

i.k a,k 



B. Density matrix and time evolution 

As stated in the introduction, our focus and goal is to study the time evolution of the distribution function of flavor 
neutrinos, at the level of free field theory at this stage. The question that we posed in the introduction and address 
here is the following: consider that at some given time the gas of flavor neutrinos and antineutrinos are described 
by a quantum statistical ensemble with a Fermi-Dirac distribution function with a fixed chemical potential for each 
flavor, namely 

n {a) (k) = — : ; n {a \k) = m * , (11.20) 

with /3 = 1/T and p a the chemical potential for each flavor. 

Such an ensemble is described by a quantum statistical density matrix which is diagonal in the Fock space of flavor 
eigenstates and is given by 



p = p (e) ® p ip) (11.21) 

with the flavor density matrices 

*(a) = e -H(H a -^q a ) ( E 22) 

Hence the initial distribution functions are given by 

<«g?«g>-^->«g?ag-»C)(fc) 

(4 a >t0<<O) = Tvp^ f3i a)i pi a) = n< a >(fc) (11.23) 

In the expressions above we have assumed that the distribution of flavor neutrinos are spin independent, of course 
a spin dependence of the distribution function can be incorporated in the description. 

Although we have stated the problem in terms of a gas flavor neutrinos in thermal equilibrium with Fermi-Dirac 
distributions, this restriction can be relaxed to arbitrary non-equilibrium single particle distributions consistent with 
Fermi-Dirac statistics. Regardless of the initial distributions the ensuing time evolution with the full Hamiltonian 
with mixing will be out of equilibrium. 



C. Cold degenerate case: 



The case of a cold, degenerate gas of neutrinos is described by the zero temperature limit but fixed chemical potential 
of the density matrix (|II.21|) with l|II.22|) . In this limit the individual flavor neutrino gases form Fermi seas "filled up" 

to the Fermi momentum kp . Consider the case of a positive chemical potential corresponding to a degenerate gas 
of neutrinos without antineutrinos at zero temperature, the degenerate ground state is given by 

\FS >= \FS > (e) <8|FS > (Al) (11.24) 

with 



\FS> ia) = Yl 



(a)t a) 

a- a- 

fc,t k,l 



|0> 



(a) 



(11.25) 
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with the flavor vacuum state |0 >( a ) annihilated by the destruction operators a ^"y fl^y The initial density matrix 
in this case is that of a pure state 



p = \FS >< FS\ (11.26) 
the distribution function of flavor neutrinos is given by 

n^(k) = 6(4 Q) - k) , n^{k) = (11.27) 
and the chemical potential is p a = u^kp). The Fermi momentum is as usual given by 



ft, p 



1 /3 

(3^)) 1/3 ^4«)(e^) = 6.1 9 ( I ^-,) ' (11.28) 



with J\f( a ^ the neutrino density for each flavor. Although the zero temperature limit is described by a pure state, this 
state is a truly many body state 

An important many body aspect of the situation under consideration can be gleaned by studying how the creation 
and annihilation operators of mass eigenstates act on the state \FS >. Consider for example the action of the 

annihilation operator a9\ on the state, to understand this question we must first obtain in terms of the creation 

Kf A /C, A 



and annihilation operators of flavor eigenstates. From equation (|II.10(l and the relation between fields given by l|II.2(l 
we find 



Cv e {k) - Sv^k) (11.29) 



and the expansion for the flavor fields given by i|II.7|) clearly indicates that if k < k F < k F , for example, then 
a^^\FS > is a superposition of states with an electron neutrino "hole", an electron antineutrino, a muon neutrino 

"hole" and a muon antineutrino. The antiparticle components of the wave function a^^\FS > is a result of the 
non-vanishing overlap between the positive energy spinors for mass eigenstates and the negative energy spinors for 
flavor eigenstates 0- 

D. Time evolution 

Within the framework of free field theory of mixed neutrinos, the time evolution is completely determined by the 
total Hamiltonian H given by eqn. I|II.14|) . 

In the Schroedinger picture the density matrix evolves in time with the full Hamiltonian as follows 

p(t) = e - iHt p(oy m (ii.3o) 

Since the full Hamiltonian H does not commute with H e , because of the flavor mixing, the density matrix does 
not commute with the Hamiltonian and therefore evolves in time. This is the statement that the initial density matrix 
(|11.21f> describes an ensemble out of equilibrium when flavor neutrinos are mixed. 

Our goal is to obtain the time evolution of the distribution functions for flavor neutrinos and antineutrinos, namely 

n^(k, t) = Trp( Q )(i) af^afl = Trp^O) ai a) x \t)ai a) x (t) (11.31) 

and similarly for the antineutrino distribution function. The initial distribution functions n^ a \k, 0) = n^ a \k) (and 
similarly for antineutrinos) given by equations l|II.23() or III. 201) for the case of an initial thermal distribution. 

It is more convenient to describe the time evolution in the Heisenberg picture wherein the density matrix does not 
depend on time and the Heisenberg field operators carry the time dependence as made explicit in eqn. I|II.36(I . 
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The free fields associated with the mass eigenstates %pi evolve in time with the usual time dependent phases multi- 
plying the creation and annihilation operators, namely 

^(fc, t ) = e lHt Mk, 0)e- lHt = Y (af e~ m ^ + e lE ^ G% ) (11.32) 

A 

The time evolution of the fields associated with flavor eigenstates, namely v a is not so simple: 



A 

where the time dependent operators ct^ (i) can be obtained by writing the flavor fields in terms of the mass 

\z . A — /c , A 

eigenstate fields using eqn. Ijll.2l) and projecting out the components using the orthogonality property given by eqn. 
(|fl.9|l . leading for example to 



ai e) (t) 

fc,A v ; 



fc.A 



/3 (e > f (t) = 

' \ y ' -fc,A 



CMk, t) + sMk, t) 



(11.34) 



The expression (|II.34|I reveals several subtle aspects which are highlighted by considering in detail for example the 
time evolution of the operator that creates electron neutrinos (a similar analysis holds for the muon neutrinos and 
their respective antiparticles) 



fe,A w L^i \ \ fc.A' fc.A' fe,A fc.A' fc.A' fc.A/ 

A' 

(Cb {1 \ , e~* E ^ G {1) ^ U<?\ +Sb (2 l x G (2 >\ U^) } (11.35) 

\ — fe,A' — fc,A' fc.A — fc,A' — fc,A' fc,A/ ) x ' 

It is a simple and straightforward exercise using the completeness and orthogonality of the respective spinor wave- 
functions, to show that the creation and annihilation operators of flavor states indeed fulfill the canonical anticom- 
mutation relations. A Fock representation of flavor states is therefore consistent and moreover needed to describe a 
quantum statistical ensemble of flavor neutrinos. 

The first line in the above expression shows that the annihilation operator for electron corresponds to the expected 
combination of creation operators for mass eigenstates multiplied by the cosine and sine of the mixing angle, but 
also multiplied by the overlap of the different spinor wavefunctions. Furthermore, the electron creation operator also 
involves the annihilation of antiparticles associated with the mass eigenstates, a feature recognized in ref . |40| . There 
are two important consequences of the exact relation (|II.35I) : 

• The amplitude for creating a mass eigenstate out of the vacuum of mass eigenstates by an electron neutrino 
creation operator is not only given by the cosine or sine (respectively) of the mixing angle, but also by the 
overlap of the spinor wave functions F^Wlli ^ . 

fc.A' fc,A 

• The electron neutrino creation operator destroys antiparticle mass eigenstates. While this aspect is not relevant 
when the electron neutrino creation operator acts on the vacuum of mass eigenstates, it becomes relevant in a 
medium where both particles and antiparticles states are populated. 

These aspects, which were also highlighted in references [H EJ E3 will be at the heart of the subtle many body 
aspects of neutrino mixing which contribute to the time evolution of the distribution functions studied below. 
The time dependent distribution functions are obtained by taking the trace with the initial density matrix 



n {a) {k, t) = Trp^^O) afUt) afcht) (11.36) 

K , A K , A 

and similarly for the other distribution functions. One can use the expression (|II.35ll for the time evolution of the 
Heisenberg field operator (and the equivalent for the hermitian conjugate), however in order to compute the time 
evolved distribution function we would need to compute the expectation value of bilinears of the field operators in 
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the flavor diagonal density matrix p(0). To do this we would have to re- write the creation and annihilation operators 
a^ ; b~ } ; etc. in the expression (III. 3511 back in terms of the creation and annihilation operators a^: /3-°^; etc.. This 

k,\ k,\ fe,A k,\ 

is obviously a rather cumbersome method. A more systematic manner to carry out this program is presented below. 
Using the expressions (|II.17III.18|) we find the following identities 

-{v a {k, t) 7° v a (k, t)) = n {a \k, t) - n (a) (fe, t) + 1 (11.37) 
' {O a (k,t){ff-k + m a )v a (k,t)) = n^(M)+n (a) (fc,t)-l (11.38) 



2u) a (k) 



Thus the computation of the distribution functions or combinations of them requires to find general expressions of 
the form 

< v e {k,t)Ov e {k,t) >=O fg < [ve{k,t)] f [v e (]e,t)] g > . (11.39) 

where the Dirac indices /, g are summed over and the averages are in the flavor diagonal density matrix ljII.21III.22p . 

Since the time evolution of the fields ipi is that of usual free Dirac field in terms of positive and negative frequency 
components, we write 

^(Ml^Aflfcle^'+A^We'^'J^tM). (11.40) 

Where we have introduced the positive and negative frequency projector operators A + (/c) and A_(fc) respectively 
which are given by 

^ -E'en? <"■«» 
AWfl-E^X-^C^w*) (,, - 42) 

ho = 7°£i(*)-7-* (11.43) 
These projection operators have the following properties, 

A^ )f (fc) = A$?(jfe) ; A^! }t (fc) = A<?(jfe), (11.44) 
AW(fc)A W (fc) = ; A { i\k)A ( i\k) = 0, (11.45) 
A^(jfc) + A w (fc) = 1. (11.46) 
We can now write the time evolution of the flavor fields in a rather simple manner by using the relations between 



the fields given by (|II.2jl and the inverse relation which allows to write ipi(k,0) in ljll.40jl back in terms of v a (k, 0). 
We find 

Mk,t) = ^Fx(k,t)[Cv e (k,0)- S Vll {k,0)] (11.47) 

Mk,t) = [CP e {k,0) - 5P M (^,0)]Fi(^t)7° (11.48) 

Mk,t) = j°F 2 (k,t)[C^(k,0) + Sv e {k,0)] (11.49) 

Mk,t) = [CP„(k,0) + S9 e {k,0)}F 2 (k,t)-/° (11.50) 

Where we have introduced the following time evolution kernels 

Fj(k,t) = 7°[A ( j 7 ' ) (fc)e-^ (fc)i + k U \ky E i^% (11.51) 

F 3 {k,t) = F 3 {k-t)f ■ 3 = 1,2. (11.52) 
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After straightforward algebra using the mixing transformation pi. 2(1 and equations l|II.47III.50|) we find the following 
result for the time evolution of the flavor fields 



v e {k,t) = T ee (fc, i)^ e (fc, 0) + T efJ- (k, t)v^(k, 0), 

v e (k,t) = v e (k,0)f ee (k,t) + P M (fc,0)T eM (A;,t), 

v,i(k,t) = T^(k,t)Vfj,(k,0) +Tfj, e (k,t)v e (k,0), 

v^{k,t) = 9^(k, 0)T w (fc,t) + H e (k,0)f^ e (k,t), 

where the time evolution operators are given by 



(11.53) 
(11.54) 
(11.55) 
(11.56) 



Tee(M) = 7° [C 2 F 1 {k,t) + S 2 F 2 (k,t)\ , (11.57) 

T^(k,t) = 7° [c 2 F 2 {k,t)+S 2 F 1 (k,t)] , (11.58) 

T eM (fc,t) = T Me = CS 7 ° |> 2 (M) - (H.59) 

f a p(k,t) = 7°J 1 ^(fc,-i)7°, (H.60) 
Furthermore since the initial density matrix is flavor diagonal, we find the following expectation values 



< [v (M)]/MM)] fl > - < [/> e (fc,o)] r K(fc,o)] s > [T ee (k,t)] rf [T ee (k,t)] ga 

+ < [^(fc ! 0)] r [z/ M (fe,0)] s > [T e/i (£,t)] r/ [T e/1 (£,t)] fl ,, 



(11.61) 



and similarly for the muon neutrino fields, where < • ■ ■ > stands for the trace with the initial density matrix. 

A noteworthy feature of the above exact expressions is that the time evolution of the flavor neutrino fields mix 
positive and negative frequency components of the mass eigenstates. Namely a flavor neutrino state is a linear 
combination of particles and antiparticles of mass eigenstates. Thus a wave packet of flavor neutrinos will necessarily 
mix positive and negative frequencies of mass eigenstates. This mixing between particles and antiparticles is a 
consequence of the fact that a flavor eigenstate is a squeezed state of mass eigenstates and viceversa[40j. 

A simple calculation yields the following expectation values in the initial density matrix 



< K(/c,0)] r [j/a(fc,0)] s > = 



E < 4>K > [^gwO + E < > 

. A A 



7 c V*(fc)-7-£ 



fia TTla o 

2u> a (k) 



[N a (k)] s 



(11.62) 

(11.63) 
(11.64) 



where n a {k); n a (k) are given by the expressions l|II.20|) and there are no flavor off-diagonal matrix elements at t = 
because the initial density matrix is flavor diagonal. 

Combining all the above results, we find the final compact form for the time dependent expectation values in eqn. 
(|II.39jl . namely 



<V e {k 1 t)Ov e {k,t)>=Tr N e (k)f ee (k,t)OT ee (k,t) +Tr ^(^(MJO^ftt) 



(11.65) 



1. Exact time evolution of distribution functions 
The exact time evolution (in free field theory) of flavor neutrinos is given by 



n (e Hk,t) = &\k,t) + J (e) (fc,i), 



(11.66) 
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where I&(k,t) and J^(k,t) are given by 



&\k,t) = 



J {e) {k,t) = 



4w e (fc) 



4w e (fc) 



Tr 



Tr 



N e (k)f ee (k, i)7°(^e + m e ) 7 °T ee (fc, t) 
N^{k)f^{k, i) 7 °(# e + m e )^°T efl {k, t) 



(M) = 1- &\k,t)- J (e) (M), 



where I^(k,t) and J^(k,t) are given by 

7W(M) = 
j (e) (M) = 



l 



4w e (fc) 
1 



Auo e {k) 

For the muon neutrinos and antincutrinos 



Tr 
Tr 



N e {k)f ee {k,t){^ e - m e )T ee {k,t) 
JV M (fc)T e/1 (M)(& - m e )T e ^{k,t) 



n^\k,t) = I^\k,t) + jW(jfe,t) 
where I^\k,t) and J^\k,t) are given by 



J<">(M) = 



l 



4w M (fc) 
1 



Tr 
Tr 



iV e (fc)f^ e (ft, t) 7 ° ( ^ + m M ) 7 °T Aie (fc, t) 



(11.67) 
(11.68) 

(11.69) 

(11.70) 
(11.71) 

(11.72) 

(11.73) 
(11.74) 



n(">(fc,t) =l-I M (k,t)-J^(k,t), 
where f»(M) and J^\k,t) are given by 



(11.75) 



4av(fc) 



4w M (fc) 



Tr 



Tr 



N^(k)f^(k, t){^ - m^)T^{k,t) 
N e (k)f^(k, *)(/£„ - m M )T Me (fc, t) 



(11.76) 
(11.77) 



The calculation of the traces is simplified by the observation that all of the different terms that enter in the trace, 
such as N a (k);T a:a >(k, t)7°; 7°T QiQ /(fc, i) can be written in the form 



7 °,4o(fc, t) - 7 • I(fc, i) + B(fc, t) =4(fc, t) + B(fc, i) 



(11.78) 



where the coefficient functions Ao(k, t): A(k, t): B(k, t) can be read off each individual term. Thus the traces in the 
terms above can be calculated by using the standard formulae for the traces of two and four Dirac matrices. 



E. Fast and slow time scales 



While the exact compact expressions above describe the full time evolution and provide a set of closed form 
expressions, they hide the fact that there two widely different time scales. These different time scales can be revealed 
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by unravelling the different contributions to the distribution functions as follows. Consider the expectation value on 
the right hand side of eqn. i|11.39|) for the case of the electron neutrino 

< [v e {k,t)] } [v e {k,t)} a > = c 2 < [&(£,t)]/ty»i (£,*)]« > +s 2 < fe(M)]/ftMM)] 9 > 

+ cs < [Mk,t)] f [Mk,t)} g + > (H.79) 

the case of the muon neutrino can be obtained from the expression above by replacing S — ► C; C — * —S. 

By writing each one of the fields ipi m terms of the positive and negative frequency contributions which evolve 
in time with the phases e^-^w* respectively, it is clear that in the products 4>i(k,t) ipi(k,t) there is a contribution 
that does not depend on time and terms that oscillate in time with the phases e T2«£;(fc)t These oscillatory terms 
arise from the interference between particles and antiparticles akin to zitterbewegung in principle do not vanish when 
the density matrix is diagonal in the flavor basis. In the general expectation values in eqn (|II.39(I these oscillatory 
terms will multiply the matrix elements of the form F- OG - , thus if these matrix elements do not vanish, these 

/c, A — k , A 

oscillatory terms are present. In the second line in eqn. l|II.79jl a similar argument shows that there are two types of 
oscillatory terms, e Ti( E i(k)+E 2 (k))t an( j e ^i(E 1 {k)-E 2 {k))t^ rpj^ f ormer ar j se f rom the interference between the particle 

and antiparticle states of different masses, while the latter from interference between particle states of different masses 
(or antiparticle). 

The combined analysis from solar neutrinos and KamLAND|45| suggest that for two flavor mixing Mf — M| = 
AM 2 ~ 7 x 10 -5 (e^) 2 and cosmological constraints from WMAP 46] suggest that the average mass of neutrinos is 
M < 0.23eV. Therefore even in the non-relativistic limit with k < Mi the ratio \E l (k)-E 2 (k)\/(Ei(k)+E 2 (k)) < 1CT 4 
and certainly much smaller in the relativistic limit k Mi . Hence because of the near degeneracy, or in the relativistic 
limit for any value of the masses, there are two widely different time scales of evolution for the flavor distribution 
functions. The longest one corresponding to the interference between particle states (or antiparticle states) of different 
masses while the shortest one corresponds to the interference between particle and antiparticle states of equal or 
different masses. This point will be revisited below. 

The evolution of the flavor (lepton) asymmetry highlights these time scales clearly and is studied below. 

III. DEGENERATE GAS OF NEUTRINOS: EVOLUTION OF FLAVOR ASYMMETRY 

The results obtained above are general and valid for any temperature and chemical potential (density). In this 
section we focus on understanding the time evolution of the flavor asymmetry n^ a '(k, t) — fS a ^(k, t) in the case of a 
cold, degenerate gas of flavor neutrinos. From equations (|11.37|) and (|11.79() we find 

n^(k,t)-n^(k,t) = ^- < V4(M)V>1(M) > +y < 4(k,t)Mk,t) > -1 

CS — -> 

+ ~y <4(k,t)Mk,t)+4(k,t)Mk,t)> (111.1) 

nto(k,t)-nV>{k,t) = y < v4(M)V>i(M) >+^r <4(k,t)Mk,t) > -1 

CS 

- — <4(k,t)Mk,t)+4(k,t)Mk,t)> (IIL2) 

The first line of the expressions above is time independent because the overlap between positive and negative 
frequency components vanishes, and the time dependence arises solely from the interference between different mass 
eigenstates. The time dependent terms (second lines in the above expressions) are opposite for the two flavors realizing 
the fact that the total charge of mass eigenstates equals that of flavor eigenstates and is time independent (see eqn. 

arm ). 

Furthermore the expectation values < ip\(k,i) ipi(k,t) > (no sum on i) are time independent (in the case of free 
field theory under consideration) since the interference term between positive and negative frequency spinors vanishes. 
The time dependence is completely encoded in the contribution that mixes the mass eigenstates. 

Therefore the time dependence of the flavor asymmetry is completely determined by the quantity 



CS / \ 
x (k,t) = — < i^l(k,t)Mk,t)+4(k,t)Mk,t)j > 



(III.3) 
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Using equations (jll.47|l - (|II.50fl . it follows that 



x(M) = ——Tr 



[N e (k) - ^(fc)][Fi(fc,t)F 2 (fc,t) + F 2 (k,t)F 1 (k,t)] 



The computation of the traces is simplified by writing 



Fj(k,t) 
Mj(t) 



co$(Ej(k)t), 
''' ahx(Ej(k)t), 



Ej(k) 



E~{k) 



sm.{Ej{k)i). 



and similarly we write 



N a (k) 

Ql 

Qa 



Q a + M a 

l 



,(«) 



(k) + 1 



k 



2u a (k) 
2to a {k) 



n^(k) - I + n^(k) 
'(k) - 1 + fjW(fe) 



For further convenience, we define 



in terms of which we obtain 



? e -#„, ; AM = M e — M M , 



(in.4) 

(III.5) 
(III.6) 

(III.7) 

(m.8) 



(III.9) 
(111.10) 

(III.ll) 
(111.12) 
(111.13) 



(III. 14) 



x(M) = — — Tr 



AM)(A(-t) + M 1 (-t)) 7 °(/ , 2 (t) + M 2 (t)) 



+ (A$ + AM)(AH) + M 2 (-i))7°(A(*) + Mi(t)) 
After some lengthy but straigthforward algebra we find 



(111.15) 



x(M) = x(fc,0)-2C 2 5 2 [(n^(fc)-n< e )(fc)) - (n^(fc)-n^(fc)) 



fc 2 + Ai!Af 2 ^ . 2 / E 1 (k)^E 2 (k) / 

1 + s 1 (fc)£; 2 (fc)J S111 1 2 



where x(&j 0) is given by 



E 1 (k)E 2 (k) 



(111.16) 



X(fc,0) = 2C* 2 S 2 [(n (e) (fc) - n^(fc)) - (n W (*0 - 



(111.17) 



The expression (|III.16|) for the time dependence of the flavor asymmetry clearly shows that neutrino mixing results 
in a time evolution of the flavor asymmetry unless the flavor asymmetry for both flavors is the same. This is obviously 
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a consequence of Pauli blocking: if the neutrino states are occupied up to the same momentum electron neutrinos 
cannot transform into an (occupied) muon neutrino state and viceversa. 

In the case of a cold, degenerate gas of flavor neutrinos (we assume here both chemical potentials to be positive) is 
given by 



n^(k) -> 9(4 Q) " k) ; n (Q) (fc) (111.18) 

If the chemical potential is different for the different flavors, the expression above shows that each wavevector mode 
will evolve with a different frequency and as a consequence of free field evolution there is no mode mixing. The 
important question is what is the time evolution of the total charge which is the integral of the flavor asymmetry over 
all momenta. This time evolution will be a result of the dephasing through the oscillations between different modes 
that are not Pauli blocked. 

We now proceed to study analytically and numerically the time evolution of the flavor charge densities q a /V with 
q a given by eqn. I|II.18|) and V the volume. We begin by defining 

u _ M x + M 2 _ AM 2 sM 2_ M | ; (nI19) 

so that Mi and M2 can be written in terms of M and AM 2 as 

*=«( 1+ IS) ■ <»-) 

We take the following as representative values for the two flavor case[iJE3 M ~ 0.25 eV and AM 2 ~ 7 x 10~ 5 (ey) 2 . 
In what follows we assume that kp > k F and introduce dimensionless variables by taking kp as the common scale, 
the opposite limit for the Fermi momenta can be obtained simply from the results below. Hence we define 

k k^ 

q = ~pr ! 9r = t~- ; T = k e F t, (111.21) 

Kp Kp 

m = — ; ,5m 2 = m\ - m 2 = 2 111.22 

k F [kpy 

(, AM 2 \ / AM 2 \ /TTTOO , 

m t = m l + w ;m 2 =fn l-_ ; (111.23) 



£1 = \Jq 2 + m{ ; e 2 = \l q 2 + m|. (111.24) 

Hence, in terms of Af (a) = (fc^) 3 /37r 2 (see eqn. l|II.28|l V we find that the time evolution of the flavor charge densities 
are given by 



qe{t) - - 6C 2 S 2 M^ (l f (r) + I s (t)) (111.25) 
tfW + §C 2 S 2 N {e) (l f {r) + I»(r)) (111.26) 



V 

qjf) 
v 



where 



If(r) 
Is(t) 



dqq 2 1- 



7711TO2 



dqq 2 1 



+ TOim 2 
£1^2 



sin 



gl +£2 

2 

£1 - £2 



(III. 27) 
(III. 28) 



We have separated the contributions from the fast (If(r)) and slow (7 s (r)) time scales as discussed in section 
(|II Efl above. In particular, as discussed above the term that oscillates with the sum £1 + £2 is a consequence of the 
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overlap between particles and antiparticles. The pre-factors that multiply the sine functions in equations ljIII.27IIII.28p 
arise from the overlap between particle- antiparticle spinors in 1111.27(1 and particle-particle, anti-particle-anti-particle 
spinors in ((111.28(1 . The overlap between particle and antiparticle spinors is non- vanishing for different masses. Similar 
contributions from the overlap between particle and antiparticle states of different masses have been found in the 
studies of refs.pollillli^. 

Since the mass eigenstates are almost degenerate or alternatively for any values of the masses in the relativistic 
limit we find 




q 2 + m\m2 
£i£2 

e = y/q 2 + fh 2 (111.29) 

with ~ 3 x 10~ 4 . Therefore the coefficient that results from the overlap between the particle and antiparticle 

spinors of different mass is given by 

q2 + mim2 o(^*)' (111.30) 



£i£ 2 V 4Af 2 E(k) ) 



and the coefficient that results from the overlap between particle-particle or anti-particle-anti-particle of different 
masses is 



2 



£l£ 2 V 4M2 E ( k ) ) 

where E(k) is an energy scale. 

Therefore the coefficient of the oscillatory term in //(t) is a factor at least of order ) ~ 1CP 7 smaller than 

that of I s (t) . Furthermore it is clear that the interference terms between particle and antiparticle average out on a time 
scale tf < 1/M whereas the particle-particle contributions evolve on a much slower time scale t s ~ M/AM 2 ^> tf. 

However, despite the fact that the coefficients of the oscillatory terms in I s {t) and //(r) differ by several orders 
of magnitude, the fact that the time evolution of I s {t) is much slower allows for a time scale within which both 
contributions are comparable. This can be gleaned from the following argument. 

The integrals for I s (t) and If(r) are dominated by the region q ~ 1. Consider an intermediate time scale so that 
the argument of the oscillatory function in If (r) is of order one, but the argument of the oscillatory function in I s (r) 

is <C 1. The contribution to the integral in i/(r) is of order fh 2 ^ j while the contribution to the integral I s (t) 

is of order 2(5m 2 T 2 ). Therefore, it is clear that even when the prefactor of its oscillatory term is small, the integrand 
of If (r) will be larger than that of I s (r) in the time domain during which 



^ > (5m 2 r) 2 =S> t < l/fh (111.32) 
In the opposite limit, for r >> l/fh the dynamics is completely dominated by I s {t). 

Fig. below displays the early time evolution of I s (t) and If(r) for < r < l/fh. It is clear from this figure that 
J/(r) averages out to its asymptotic value on a short time scale r~ 1 (t ~ 1/kp) and that I s (t) begins to dominate 
the dynamics on time scales r > l/fh as discussed above. In the case of Fig.Q, with k% » M the time scale of 
averaging is t ~ 1/kp, but for kF « M it would be of order 1/M. 

In terms of dimensionful quantities the inequality obtained in eqn. 1(111.32(1 above translates into t < 1/M. With the 
current estimate M ~ 0.25 eV the analysis above suggests that the particle-antiparticle interference is dynamically 
relevant during time scales t < 10~ 15 s although this time scale is comparable to the expansion time scale at the time 
of the electroweak phase transition, it is far shorter than the time scales relevant either for primordial nucleosynthesis 
or for dynamical processes during the collapse of supernovae or neutron star cooling. 

While the behavior of I s {r) and If(r) as a function of r must in general be studied numerically, the long time limit 
can be extracted analytically. 

The asymptotic long time behavior of I s (t) and If(r) is determined by the end points of their integrands, in 
particular for momenta near the Fermi surface. Two relevant cases yield the following results 
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FIG. 1: J s (r) and I s {r) for fcf. = 100 eV ; fc£ = ; M = 0.25 eV ; AM 2 ~ 7 x 10~ 5 (eV) 2 vs. r. For these values l/fh = 400. 



Relativistic Limit: max(k F ,k F ) ^> Mi,M 2 



Is(r) 
I fir) 



dqq 1 



£l£ 2 



Sm 2 i 



sin 



(577J 2 



(jf 2 + Tfl\Ul2 \ 1 



— (mi - ma)' [sin(2r) - sin 2g r r] + 



1 



where <5m 2 is defined by equation l|III.21|) along with the other dimensionless variables. 
• Non-Relativistic limit: k c F ,k F < M 1: M 2 



Is{r) 



1 / „2 

2 



£l£ 2 



m\m,2 



(mi - m 2 )(l 



g r sm 

Is(r) = I [ dq,f ( 1 + 



(mi — m 2 ) t 
1 



(mi - m 2 )(l 



2mim 2 



)r 



2miTO 2 ' 



O 



q z + mim 2 
£l£ 2 



[1 — cos[(toi — m 2 )i 



(mi - m 2 )r 
for ^> 1 

TOlTO 2 

for (mi-m 2 )r <<; x 
mim 2 



(111.34) 



(111.35) 



I fir) 



I fir) 



dqq" 1 



q + mim 2 \ (mi — m 2 ) 2 



£l£ 2 



/ mim 2 (mi + m 2 ) r 



(mi + m 2 )(l 



1 



2mim 2 ' 



<7^ sin 



(mi + m 2 )(l 



2mim 2 



for (mi+W2)T »l 



mim 2 



rfgg 2 1- 



g 2 + mim 2 

£l£ 2 



[1 — cos[(mi + m 2 )r]] ; for 



(mi + m 2 )n 
mim 2 



< 1 



(111.36) 



In both cases, the flavor asymmetry density at asymptotically long time is given by 



-(q e (t) - q^(t)) 



cos 2 (26>) +G(l/t) 



(111.37) 



The power law fall-off is a consequence of dephasing between different flavor modes that are not Pauli blocked. 
Fig.© displays the slow contribution I s (t) and its asymptotic limit given by eqn. (|III.33(I hi the relativistic case. 



17 




0.5 

a 
T 

"5f o.o - 



m 



20 40 



80 100 



FIG. 2: 7 s (r) and [J s (r) - J 3 (oo)] x (fir) for fcj. = 100 eV ; fe£ = ; M = 0.25 el/ ; AM 2 ~ 7 x 10~ 5 (eV) 2 vs. fir, with 

n = <W = 7 x lo -9 . 



IV. DISTRIBUTION FUNCTIONS OF NEUTRINOS AND ANTINEUTRINOS 



The distribution functions are given by equations (|II.66|) - (|II.77|I for which after lengthy but straightforward algebra 
we find the following expressions 



I (e) (k,t) = n^(k)-2n {e) (k)A(k,t)- \n^(k) - (l - n^(kj)] B(k,t) (IV. 1) 

ujiik) L V /J 



J {e) {k,t) 



C 2 S 2 (M 1 - A/ 2 ) 2 + MiM 2 + k 2 
u e {k)w tJl {k) 



lM(k) - [1 - n ( ^(fc)]) + (n^\k) + [1 - nW(jfe)]) 



nW(fc)- (l-nM(fc)) 



C(M) 



where n^ e \k) and fi^ e ^(fc) are the initial distribution functions given by eqn. flll.20f) and 



A(k,t) 

(IV.2) 



A(k,t) = C 2 S 2 



_ k 2 +M 1 M 2 \ . 2 ( E 1 {k)+E 2 {k) 
E 1 (k)E 2 (k) |sm 1 ' ' 1 



k 2 + M 1 M 2 \ a .^ 2 (Ei{k) - E 2 {k) 
E\(k)E 2 {k) 



sm 



B(k,t) = C*S*(Mi - M 2 ) 



C(M) = C 2 5 2 (Mi-M 2 ) 2 
2C 2 S 2 - 1 



Ei (AO 
C 2 S 2 



sin(£i(fc)t) -—-sm(E 2 (k)t) 
E 2 (k) 



sin 2 (S 1 (fc)i) + — — sin 2 {E 2 (k)t) 



E 1 {k)E 2 {k) 



E 2 {k) ^ " E 2 (k) 
sin(E 1 (k)t)8Wi(E2(k)t) . 



(IV.3) 
(IV.4) 
(IV.5) 

(IV.6) 



The expressions for I^(k,t) and J&(k,t) are obtained from those for 1^ (k, t) and (k, t) above by the replace- 
ment 



n {e \k) < — ► [1 - n w (fe)] ; < — ► [1 - fi ( ^(fc)] 



(IV.7) 



Finally the expressions for (fc, i); (fc, i); j00 (fc, t); (k, t) are obtained from those for the electron neutrino 
by the replacement 



n {e \k) < — ► n^(fc) ; fi (e) (fc) < — ► n M (k) ; w e (A) < — ► w M (Jfc) ; C 2 



S 2 



(IV.8) 
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These dynamical factors A(k, t);B(k,t);C(k,t) are determined by the time evolution while their pre-factors in 
the expressions for the distribution functions are determined by the initial state. The dynamical factors clearly 
reveal again the different time scales. Terms that feature the contributions e^ 2 *^ 1 - 2 )*; e ±l i E ^+ E ^) t arise from particle- 
antiparticle interference and their contribution is proportional to (AM 2 /M 2 ) and those that feature e ± J ( s i- B 2)t ar i se 
from particle-particle (or anti-particle- anti-particle) interference. We can find the asymptotic distribution functions 
at long time by averaging the oscillatory terms over a time scale longer than the longest scale ~ M / AM 2 . This time 
averaging procedure leads to 



A(k,t) 



C 2 S 2 



B(k,t) = C(k,t) = -C*S*(Mi - M 2 y 



E 2 (k) E 2 {k) 



(IV.9) 
(IV.10) 



The above expressions are exact and therefore valid for any value of the neutrino masses M 1; M 2 . However, the 
most recent compilation ^3 of data suggests that in the two flavor case the mass eigenstates are almost degenerate 
with AM 2 ~7x 10 -5 (eV) 2 and the most recent cosmological constraint from WMAP|46| suggests that the average 
value of the mass M is < 0.25 (eV). In terms of the M and AM 2 introduced in eqn. (jIII.19|l . we find 



m P = M 



AM 2 

1 + —^r C0S(2#) 

AM 2 K ' 



; m M = M 



AM 2 

1 =— cos(26») 

AM 2 y ' 



(IV. 11) 



In terms of the small ratio AM 2 /M 2 ~ 10 3 wc find the average of the distribution functions over the longest time 
scale to be given by 





= n (e \k)- 


2C 2 S' 




= n (e) (fc)- 


2C 2 S : 




= n ( ^(fc)4 


2C 2 S 




= n^(fc) + 


2C 2 S 



,(/' 



} (fc)) - n[k,n {a \n^] 
--R[k in ^ a \n {a ^] 

-n[kM a) ,fi {a) } 

-TZ[k,n (a) ,n^] 



(IV.12) 

(rv.13) 

(IV. 14) 
(IV. 15) 



with 



n[k,n (a \n ( ^] 



k 2 M 2 /AM 2 " 
E 4 (k) V AM 2 

( / AM 2 * 



C 2 S 2 



AC 2 S 2 (n^(k)+n {e \k) - 1) + (n w (fc) + n M (k) - 1) 



-O 



V AM 2 



E{k) = Vk 2 + M 2 



(IV.16) 



(IV. 17) 



The term TZ[k, n^ a \ n^} arises from the overlap between particle and antiparticle spinors which features the small 
quantity (fj^). 

Flavor pair production and normal ordering: 

The expressions (|IV.12IIV.15ll with that for the corrections given by eqn. (|IV.16(1 point out an important and 
subtle aspect of the dynamics of mixing. Consider that the initial density matrix is the flavor vacuum, namely set 



.(e) 



(ft) 



i(e) 



(k) 



(k) = 0. The asymptotic limit of the distribution functions ljIV.12IIV.15p is given 



to lowest non trivial order in the ratio AM / M by 



(k, oo) 



(fc, 00 ) 



,00 



(k, oo) 



5 00 



(k, oo) 



k 2 M 2 



( 



AE 4 (k) \ AM 2 
< ( AM 2 \^ 



sin 2 (2(9)(l + sin 2 (2(9)) 



O 



V AM 2 



(IV. 18) 
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This result clearly indicates that the time evolution results in the creation of particle-antiparticle pairs of electron 
and muon neutrinos. This is of course a consequence of the non- vanishing overlap between positive and negative 
energy spinors which results in that a destruction operator for flavor neutrinos develops a component corresponding 
to a creation operator of antineutrinos during time evolution, and viceversa. In leading order in the degeneracy, the 
typical momentum of the pair created is k ~ M therefore these are typically low momentum pairs of flavor neutrinos. 

Furthermore a remarkable aspect of this pair production process via neutrino mixing is that the distribution 
function of the produced particles falls off very slowly at high energies, namely n proc i(k, oo) cx 1/fc 2 . As a result 
there is a divergent number of pairs produced as a consequence of mixing and time evolution. Since the particles 
and antiparticles are produced in pairs, the flavor charge vanishes, but the individual distribution functions feature 
a contribution from the pair production process. A normal ordering prescription must be appended to subtract the 
infinite number of particles created, however unlike normal ordering in the usual free field theory, which subtracts a 
constant, in the case of mixing such normal ordering requires a subtraction of a distribution function. 

This is a novel and subtle phenomenon, flavor pair production which is a direct many particle consequence of mixing 
and oscillations. Since this phenomenon is a consequence of the interference between particle and antiparticle states 
is suppressed by the small quantity (AM 2 / M) 2 . 

Regardless of whether this phenomenon of flavor pair production has any bearing on the cosmology and or astro- 
physics of neutrinos, it is a genuine many body aspect inherent to the field theory of neutrino mixing that deserves 
to be studied in its own right as a fundamental aspect of the field theory of mixing. 

Off-diagonal densities: Even when the initial density matrix is diagonal in the flavor basis and therefore there are 
no off-diagonal initial correlations, these develop upon time evolution as a consequence of flavor mixing. Following 
the same steps described above for the distribution functions, we find the off-diagonal density to be given by the 
following expression 
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with A(k,t) given by eqn. ljIV.3|) . The expressions for the distribution functions and the off-diagonal density can 
be simplified by expanding the coefficients of the oscillatory functions up to leading order in the small quantity 
(AM 2 /M 2 ). We find 



n^(k,t) = n {e \k) - (n {e \k) - n^\k))2C 2 S 2 [1 - cos[(^i(fc) - E 2 (k))t]] + O (IV.20) 

\ 4M 2 / 

the other distribution functions may be found from the expression above by the replacements in eqns. (|IV.7IIV.8|I . 
Their time averages over the long time scale coincides with the leading expressions in eqns. (|IV.12IIV.15|) . The 
off-diagonal density simplifies to the following expression 
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The terms of O ( 4M2 ) again involve terms that oscillate with the sum of the frequencies corresponding to particle- 

antiparticle interference as well as terms that oscillate with the difference of the frequencies arising from the overlap of 
the particle (or antiparticle) spinor wavefunctions for different masses. The analysis that was presented for the same 
type of contributions in / s (r); If(r) above highlight that the particle- antiparticle interference becomes subdominant on 

jT I m the approximations ijlV.201 IIV.21|) determine the dynamics 

of the distribution functions and the off-diagonal correlator in leading order in the small ratio 4rp" f° r t >> l/M. 



time scales t > l/M. Hence the first terms O 



A. Equilibrated gas of mass eigenstates 

Although we have focused on the case in which the initial density matrix is diagonal in the flavor basis, for 
completeness we now study the case in which the initial density matrix describes an ensemble of mass eigenstates in 
equilibrium. Therefore this initial density matrix is diagonal in the mass basis and commutes with the Hamiltonian. 
This situation thus describes a state of equilibrium in which the occupation numbers do not evolve in time (in the 
non-interacting theory). In this case we find 
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where n^ l \k) are the occupation numbers of mass eigenstates assumed to depend only on the energy. Just as we did 
in our previous analysis it proves convenient to write the above correlator in the following form 
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2Ei{k) 

Since the density matrix commutes with the full Hamiltonian, the distribution functions of the flavor eigenstates 
do not depend on time. Following the procedure detailed above we find the following results 
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Using the relations given by eqn. (|IV.11|) we find to leading order in AM 2 / M 2 
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V. "EFFECTIVE" (FREE) FIELD THEORY DESCRIPTION: 



Let us summarize the lessons learned in the analysis of the previous section in order to establish a set of criteria 
with which to develop an effective description of the dynamics in the case in which the mass eigenstates are nearly 
degenerate as confirmed by the experimental situation or in the relativistic case. 

• For nearly degenerate mass eigenstates there is a hierarchy of scales determined by i) kp or temperature (T), 
ii) the average mass M and iii) the mass difference M\ — M-z- The experimental situation seems to confirm the 
near degeneracy with \M\ — M%\ <C M, therefore at least two scales are widely separated. Furthermore if kF 
and or T (temperature) are such that kp;T ^> M which describes a relativistic case, then all three scales are 
widely separated with the hierarchy kp,T 3> M ^S> \M\ — M%\. The dynamics studied above reveals all three 
scales. 

• The time evolution of the distribution functions, flavor asymmetry and off-diagonal correlators all feature terms 
that oscillate with the frequencies E\(k) + E2(k), 2Ei^(k), and also terms which oscillate with the difference 
Ei(k) — i?2(fc). The former arise from the interference between particle and antiparticle states of equal or 
different masses and determine the short time scales t < 1 /M , while the latter arise from interference between 
particle states (or antiparticle states) of different masses and determine the long time scales t > M/AM 2 . The 
terms that oscillate with the fast time scales average out on these fast scales and their coefficients arc of order 
AM 2 I'M 2 and hence small in the nearly degenerate case. These coefficients result from the overlap between 
positive and negative energy spinors of slightly different masses. The coefficients of the terms that oscillate on 
the long time scale are of 0(1) and result from the overlap between positive energy spinors (or between negative 
energy spinors) of different masses. 

• The contributions to the distribution functions and off-diagonal correlators from the terms with fast and slow 
oscillations are comparable within the short time scale t < 1/M but for times longer than this scale the contri- 
butions from the terms with fast oscillations are suppressed with respect to those with slow oscillations at least 

by of^C 



\ 4M 2 

We seek to obtain a description of the oscillation dynamics on scales much larger than 1/M when the contribution 
from the fast oscillations have averaged out to quantities that are proportional to powers of the small ratio ^£j 2 and 
can therefore be neglected in the nearly degenerate case. 

In the nearly degenerate case AM 2 /M 2 -C 1 the masses m e , m M , M±, M 2 ~ M (see eqns. (|IV.ll(lllIL20(l ). thus_in 
order to isolate the leading order terms as well as to understand corrections in the degeneracy parameter AM 2 /M 2 
it proves convenient to expand the positive and negative energy spinors in terms of this small parameter. A straight- 
forward computation in the standard Dirac representation of the Dirac gamma matrices leads to the following result 
for the flavor positive and negative energy spinors (see eqn. ill A 
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with 



E(k) = \Jk 2 + M 2 



(V.2) 



and the upper sign corresponds to a — e and the lower sign to a = [i. The spinors U? A > A are positive and negative 
energy solutions respectively of the Dirac equation with mass M with unit normalization. Similarly for the positive 
and negative energy spinors associated with the mass eigenstates F~ ' ; G g (see eqn. HlHOf) ). we nn d 
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with the same spinors A ; A , where the upper sign corresponds to i = 1 and the lower sign to i = 2. 

It is clear from the approximations l|V.l(l and (|V.3j) that the overlap between positive and negative energy spinors 

of different masses is O ^ ^jp • For times much larger than the fast time scale, the corrections to the spinors are 
subdominant and can be neglected and the fields associated with the flavor and mass eigenstates are expanded as 



MM) = £ ( a gJ(t) W^ + Z^t (t) V_ S , A ) +0 (^) (V.4) 

^(M) = E (««> W SiA e-^«« + 6^ V_ S , A +0 (^) . (V.5) 

We can now find the relation between the creation and annihilation operators of flavor states and those of mass 
eigenstates by using eqn. (|II.2J| . to leading order in the degeneracy parameter we find 
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where we have neglected terms of O f ^ps ) i an d similar relations hold for the annihilation operators of the respective 

antiparticles fi^ 1 ' (t) . It is clear that the approximations leading to the relations (|V.6(1 and (|V.7(I are more generally 

valid not only in the nearly degenerate case but also in the relativistic case k ^> regardless of the value of the 
mass difference, since in this case the common spinors are those of massless Dirac fermions in all cases. 

In this approximation, the evolution equation for the Heisenberg operators (t) does not follow directly from any 

Dirac equation, but can be obtained straightforwardly by taking time derivatives of these operators in eqns. I|V.6IV.7|) 
and using the relations (|V.6IV.7|I to re-write the result in terms of the operators themselves. In the leading order 
approximation particles and antiparticles do not mix since the overlap between the spinors . and V_j A vanishes 
(in free field theory) and a straightforward calculation leads to the following equations of motion 
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and a similar equation of motion for the annihilation operators for flavor antiparticles (3j? '}(t)- These equations 

of motion look to be the familiar ones for neutrino oscillations0-0, IT2I flfl I22I l23l |. but these are equations for 
the Heisenberg field operators, rather than for the single particle wave- functions. Once the time evolution of the 
operators is found, we can find the time evolution of any multiparticle state. Furthermore the regime of validity of 
these equations is more general, they are valid either in the nearly degenerate case AM 2 /M 2 <C 1 for any value of 
the momentum, or in the relativistic limit for arbitrary value of the masses provided that fc ^> M 1; M 2 . 

Inverting the relation between the operators for flavor and mass states at the initial time, namely writing the 



operators a- in terms of a- (0) using eqns. I|V.6IIV.7|) at t = 0, we find (again to leading order) 

k > X Jc » X 



al e) (t) = al e) (0) [c 2 e- iEl «* + ^-^(fc)*! + sc a M {Q) 

k.X h \ x / h \ v 



k,X' 



k,X y ' 



k,x y ' 



(j2 e -iE 2 (k)t , g2 e ~iE 1 (k)t 



SC al e ^(0) 



-iE 2 {k)t 



,-iE 2 (k)t 



-iE 1 (k)t 



-i_Ei(/c)t 



CV.11) 

(V.12) 



For the antiparticle operators we find the same equations with 

k.X 



The Heisenberg field operators given by eqns. HV.lllV.12p (and the equivalent for the antiparticle operators) are 
the solutions of the equations of motion l|V.8|) . 

The time evolution of the distribution functions in an initial density matrix that is diagonal in the flavor basis 
follows from a straightforward calculation using the above time evolution. We find 
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The distribution functions for antiparticles to leading order is obtained from the above results by the replacements 
A straightforward calculation following the above steps leads to the result 
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The results (|V. 1 3|) and (|V.15J| reproduce the leading order expressions found in the previous section, eqns. 
(|IV.20IIV.21|) . Thus this "effective" free field theory description reproduces the leading order results either in the 
nearly degenerate case AM 2 <C M 2 or in the relativistic case. Furthermore either the effective equations of motion 
l|V.8Jl or alternatively the time evolution ljV.lHV.12JI (and those for antiparticles) lead to a set of closed evolution 
equations for bilinear s. These are most conveniently written by introducing a fiducial spin S = (S x , S y , S z ) with the 
following components 
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and a fiducial magnetic field B = (B x , B y , B z ) with components 

~B(k) = 2O(fc)(0, - sin(20), cos(20)) 
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in terms of which the equations for the bilinears are akin to the Bloch equations for a spin S precessing in the 
magnetic field B namely 
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The antiparticle operators obey independently a similar set of equations. To leading order in AM 2 /M 2 there is 
no mixing between particles and antiparticles (suppressed by two powers of this small ratio) , therefore the number of 
electron plus muon neutrinos is conserved independently of that for antineutrinos, namely 



and similarly for the operators ?■ The set of equations above, for Heisenberg operators is akin to the equations 
of motion for the "single particle" density matrix obtained in ref.[2^|, which are equivalent to those investigated in 

refsmmnsmin. 

In the study of synchronized oscillations[2^|23,|^,|23, a self-consistent Hartree-Fock approximation is introduced 

which leads to a Bloch equation like l|V.20J) but where the magnetic field B acquires a correction from the self- 
consistent Hartree terms which arise from forward scattering off neutrinos in the medium. 

This effective formulation neglects the dynamics of flavor pair production discussed above since such phenomenon 
is suppressed by two powers of the small ratio AM 2 /M 2 . 



A. Propagators: non-equilibrium correlation functions 



While the set of equations of motion ljV.81) and (|V.20() are reminiscent of those for the single particle wave functions 
and the single particle density matrix, in fact there is more information in the "effective" free field theory description 
afforded by the operator equations l|V.8|) and (jV.20|l combined with the field expansion l|V.4|l . In particular, inserting 
the solution of the equations of motion (|V.11I IV.12|I (and the similar ones for the antiparticles) into the expansion 
(|V.4f) for the field operators allow us to obtain any correlation function in the free field theory at equal or different 
times. These are the building blocks of any systematic perturbative expansion of processes of weak interactions. In 
particular the Feynman propagators, which are an essential ingredient in any calculation that involves neutrinos are 
given by 
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where the expectation values are in the initial density matrix, which is taken to be diagonal in the flavor basis in the 
present discussion. 

The correlation (Wightmann) functions that enter in the Feynman propagator are found by using the leading order 
expansion ljV.4l) with the time evolution of the creation and annihilation operators given by eqns. HV.lllV.12p and 
similar ones for f3 [ *\f). With the purpose of highlighting the fast and slow time scales in the propagators, it is 
convenient to introduce the following functions that evolve on the slow time scale 
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in terms of which the Heisenberg creation and annihilation operators of flavor states are written as follows 
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and similarly for the antiparticle Heisenberg operators /3^"Ht)- 

A straighforward calculation of the Wightman functions yields the following results 
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The Wightman function for the muon neutrino is obtained from that of the electron by the replacement 
n,( e \k), n^ e \k) — > n^(k),n^(k) , and f k < — ► f k . The off-diagonal Wightman functions are given by 
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the other off-diagonal Wightmann function is obtained from the one above by replacing < — > n^' a^c? f k < — > 

We have specifically separated the "fast" evolution, encoded in the exponentials e I±B ( fc )(*~* ) and the "slow" evo- 
lution encoded in the functions f k ;g k which oscillate with the small frequency O(fc) ~ AM 2 /2E(k). We emphasize 
that the propagators above are functions not only of the difference (t — t') but also of the sum (t + 1') which reveals a 
truly non- equilibrium evolution. The manifest lack of time translational invariance reflects the fact that the density 
matrix which is diagonal in the flavor representation does not commute with the time evolution operator. 

The discussion at the beginning of this section points out that these propagators are valid on time scales t, t' » 1/M, 
for which the corrections arising from the interference between particle and antiparticle can be neglected. Therefore 
the correlation functions obtained from the effective field theory must be understood as being averaged over the fast 
time scales and their validity is restricted to slow time scales. 

The free field theory propagators obtained above provide the main ingredients to carry out a study of the weak 
interactions in a neutrino background in a loop expansion. 



VI. CONCLUSIONS AND DISCUSSIONS: 



Our focus was to study the evolution of a dense and or hot gas of flavor neutrinos as a consequence of oscillations 
and mixing. The goal was to establish an understanding of the dynamics directly from the underlying quantum field 
theory, beginning with the simplest case of free field theory and restricted to the two flavor case. 

Such study leads to a deeper understanding of the various approximations invoked in the literature as well as 
recognizing the potential corrections. Even at the level of free field theory, which must be the starting point of 
any program to study the physics of oscillations and mixing in the weak interactions, this study reveals a wealth of 
dynamical phenomena that has not been explored before within the context of neutrino oscillations in a medium with 
neutrinos at finite density and temperature. 

The most salient aspects of our study are the following: 

• A hierarchy of time scales emerges associated with different interference phenomena. Oscillations on fast time 
scales t < 1/M are associated with the interference between particles and antiparticles while oscillations on 
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slow time scales t > M/AM 2 arise from the interference between particle (or antiparticle) states with different 
masses. Observationally the situation for two flavors is that of near degeneracy, which entails that these time 
scales are widely separated. Furthermore in the relativistic limit with typical energy E 3> Mi , M2 there is an 
even shorter time scale t ~ 1/E. 

• The terms that oscillate on fast scales feature coefficients that are determined by the overlap of positive and 
negative frequency wave functions of different masses. In the relativistic limit or in the case of near degeneracy as 
suggested by the recent observations, these terms are of order (AA/ 2 / M 2 ) 2 ~ 10~ 6 (or smaller in the relativistic 
case), while the coefficients of terms that oscillate on the slow scales are of 0(1) in terms of this ratio. During 
the short time scales both contributions are comparable, but for t >> 1/M the contribution from the overlap 
between particle and antiparticle states becomes subdominant being at least a factor (AM 2 / M 2 ) 2 ~ 10~ 6 
smaller than the oscillations on the slow time scale. For the values of M consistent with the recent bounds 46] 
the scale for fast oscillations is ~ 10~ 15 s these are clearly too fast for relevant processes during BBN or neutrino 
processes in astrophysics, but may be relevant for early universe cosmology. Of course this possibility requires 
further and deeper studies. 

• An initial flavor asymmetry relaxes to equilibrium via dephasing between modes that are not Pauli blocked with 
a power law \/t on slow time scales t > fcp/Am 2 in the relativistic case Uf » M. We have obtained exact as 
well as approximate expressions for the time evolution of the distribution functions and off diagonal densities 
and discussed their asymptotic behavior, all of which display Pauli blocking between different flavors (see eqns. 
(|IV.12IIV.15|) . For completeness we have also studied the case of an equilibrated gas of mass eigenstates which 
describes a situation of equilibrium in absence of interactions. The non-equilibrium oscillation dynamics leads 
to the production of particle- antiparticle pairs of flavored neutrinos with typical momenta k ~ M . Since this 
phenomenon is a direct consequence of the overlap between particle and antiparticle states the pair yield is 
suppressed by the factor (AM 2 / M 2 ) 2 . 

• The wide separation between the different time scales allows to describe the dynamics on the longer time scales 
in terms of an "effective" theory. In this effective description the Heisenberg creation and annihilation field 
operators for flavor neutrinos and antineutrinos obey the familiar Bloch type equations and the spinor structure 
is common to both flavors as well as the mass eigenstates. This effective description allows to obtain in a simple 
manner the dynamics of the distribution functions, off diagonal correlation functions and the non- equilibrium 
propagators, all of which must be understood as an average over the fast time scales and valid only on the slow 
scales. 

While we have focused on the evolution of a gas of flavor neutrinos as an initial value problem we have not 
discussed how the initial state is "prepared" . This is an important aspect of the physics of neutrino mixing and 
the weak interactions, since weak interactions only produce flavor states the initial state (or density matrix) must 
be "prepared" by weak interaction processes that occur on time scales much shorter than those in which such state 
will relax either via collisions or by oscillations. Clearly we have nothing to say yet on this aspect which deserves a 
thorough study. 

Another aspect that deserves attention is that of the corrections to the "effective" theory described above. These 
corrections entail powers of the ratios that are small either in the nearly degenerate case or in the relativistic limit. 
In perturbation theory in the weak interactions, these "small" corrections could conceivably be comparable to pertur- 
bative corrections in Gf the Fermi coupling, in which case the terms neglected in the effective theory must be kept 
on the same footing as the contributions in the weak coupling in the perturbative expansion. Clearly such possibility 
must be evaluated for the particular situation under consideration. 

While we have focused on the dynamics in free field theory, the results will likely be valid in the interacting case 
in the case of a low density neutrino gas (or low temperatures). Under these circumstances the corrections to the 
evolution equations associated with forward scattering off the neutrino background (mean field), which is of order Gf 
would be much smaller than AM 2 /M and the free field theory results for the evolution of the asymmetry may very well 
be valid. Furthermore, the weak interactions only affect the left handed neutrinos but not the right handed neutrinos 
which will oscillate as in a free field theory. The mass term will then entangle the oscillations of the right and left 
handed components. Such a process will be suppressed in the relativistic limit but may introduce yet another scale. 
The intriguing phenomenon of flavor pair production, a many body feature intrinsic to the field theory of neutrino 
mixing and oscillations. While it is not clear to the authors whether such phenomenon could have potential bearing 
in cosmology and astrophysics, it certainly is part of the fundamental aspects of neutrino mixing and oscillations and 
deserves further study. 

We are currently studying these and other possible scenarios including interactions. 

Having understood the regime of validity of the effective "long time" theory as well as having obtained the necessary 
non-equilibrium propagators we expect to address the issue of the propagation of neutrinos in a dense and or hot 
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medium, including a neutrino background including not only forward scatterin g bu t also collisional processes P. Ifl, l2~il 
I2M l2fj| by implementing the methods of non-equilibrium quantum field theory]!?! . 
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